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ABSTRACT 


This  report  deals  with  the  development  and  utilization  of  a  maximum 
likelihood  tracking  algorithm  designed  to  handle  a  single  diffusing 
target.  The  tracker  is  required  to  accept  or  reject  each  of  a  sequence 
of  discrete  position  reports,  some  of  which  are  false  alarms.  Average 
tracking  time,  E(T),  is  defined  and  used  to  evaluate  tracker  performance. 
The  effects  of  tracker  memory  are  examined  and  a  Kalman  filter  is 
developed  to  handle  the  effects  of  measurement  error.  Simulation 
results  are  included. 
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I,  INTRODUCTION 


The  treatment  of  false  alarms  is  central  to  the  development  of  any 
useful  set  of  tracking  decision  rules.  The  purpose  of  this  report  is 
to  examine  the  effectiveness  of  a  specific  tracking  algorithm  to  be 
used  in  the  presence  of  false  alarms.  A  maximum  likelihood  decision 
rule  will  be  developed  and  used  in  conjunction  with  an  expandable 
tracker  memory.  The  primary  measure  of  tracker  effectiveness  will  be 
average  tracking  time,  E(T) . 

A.  BACKGROUND 

Every  tracking  system  must  in  some  way  account  for  the  possibility 
of  receiving  both  real  and  false  target  position  reports.  Tracking 
radars  typically  employ  a  'gate'  criterion  that  rejects  position  re¬ 
ports  which  differ  significantly  from  what  target  dynamics  might  per¬ 
mit.  The  gate  helps  to  maintain  track  in  a  relatively  noisy  environ¬ 
ment  at  the  expense  of  occasionally  rejecting  a  real  detection. 
Surveillance  systems  with  data  rates  considerably  below  those  of 
tracking  radars,  such  as  might  be  found  in  an  ASW  environment,  must 
also  deal  with  false  detections,  and  these  will  be  the  subject  of 
this  report. 

Tracking  systems  with  low  data  rates  must  necessarily  take  a  dif¬ 
ferent  approach  than  a  typical  aircraft  tracking  system  might  employ. 
Most  rapid  rate  tracking  systems  are  judged  by  their  ability  to 
recognize  and  differentiate  new  tracks,  crossing  tracks,  split  tracks, 


and  fading  tracks.  They  must  be  able  to  process  large  amounts  of  in¬ 
formation,  quickly  and  efficiently,  to  produce  precise  target  position 
estimates  on  a  real  time  basis.  In  his  article  'An  Optimal  Data  Asso¬ 
ciation  Problem  in  Surveillance  Theory*  Robert  Sittler  developed  a 
maximum  likelihood  model  for  just  this  purpose. 

In  a  situation  where  the  data  rate  is  low,  the  emphasis  shifts  from 
how  well  we  track  to  whether  we  can  track  at  all .  We  will  look  at  a 
particular  class  of  trackers  that  maintains  an  estimate  of  target  posi¬ 
tion  at  all  times,  and  for  which  it  is  therefore  possible  to  define 
tracking  error  at  all  times. 

Figure  1  illustrates  the  results  of  three  separate  tracking  se¬ 
quences  using  the  maximum  likelihood  memoryless  tracker  developed  in 
this  report.  In  the  figure,  error  is  measured  in  units  of  length  and 
used  to  represent  the  difference  between  the  target's  real  position 
and  the  tracker's  estimated  position  of  the  target  at  any  particular 
time.  Due  to  the  discrete  nature  of  the  target  position  reports  used 
in  the  simulation,  error  is  measured  accurately  only  at  the  positions 
marked  with  a  •.  In  all  cases  the  error  initially  fluctuates  up  and 
down,  but  finally  grows  very  rapidly.  We  wish  to  identify  a  "tracking 
time"  that  marks  the  end  of  the  fluctuation  region,  and  to  define  a 
tracking  scheme  that  makes  the  average  "tracking  time"  as  large  as 
possible.  As  figure  1  shows,  a  precise  definition  of  tracking  time 
can  be  illusive,  and  we  will  therefore  define  a  conservative  measure 
later  in  this  chapter. 
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□  -  First  mistaken  ID  in  track 


Time  vs  error  plot  for  3  sample  runs  of  the  maximum 
likelihood  memoryless  tracker. 


3.  PROBLEM  DYNAMICS 


We  consider  an  abstract  two-dimensional  version  of  the  tracking 
problem.  Real  target  position  reports  arrive  randomly,  according  to 
c  Poisson  process  with  data  rate  X  [units  are  (time”^)].  In  a  process 
of  this  kind,  reports  are  equally  likely  to  arrive  at  any  point  in  a 
fixed  interval.  It  is  not  possible  to  identify  real  detections  by 
their  time  of  arrival  or  the  interval  between  arrivals. 

Target  motion  is  assumed  to  be  due  to  normal  diffusion  with  no 
drift.  A  diffusion  constant  D  [units  are  (length  x  time  )]  will  be 
used  to  characterize  a  single  target's  motion.  Motion  of  this  type 
can  be  modeled  with  a  bivariate  normal  distribution.  In  a  situation 
with  no  false  alarms  or  measurement  error,  the  target  density  would 
have  the  physical  shape  of  a  tall  mountain  centered  on  the  last  target 
position  report.  The  mountain  spreads  gradually  in  all  directions  as 
the  interval  between  detections  increases,  much  like  a  mountain  of  snow 
melting  under  the  sun.  The  motion  of  Fleet  Ballistic  Missile  Sub¬ 
marines  on  patrol  is  often  modeled  with  this  type  of  distribution. 

Tracking  systems  are  often  prone  to  some  inherent  measurement 
error.  Many  systems  will  have  an  internal  bias  in  one  direction, 
around  which  a  detection  might  be  likely  to  occur.  In  this  report 
we  will  consider  the  case  in  which  that  type  of  bias  can  be  recognized 
and  accounted  for.  We  will  assume  that  the  reported  target  position 
has  a  circular  normal  measurement  error  with  standard  deviation  a 
[units  are  (length)  in  each  direction],  and  will  be  centered  on  the 
actual  target  position. 


In  most  tracking  systems  there  exists  the  possibility  of  a  detec¬ 
tion  occurring  in  a  location  where  no  target  exists.  In  a  radar  sys¬ 
tem  it  may  be  due  to  a  random  electrical  disturbance,  or  with  sonar, 
due  to  a  transient  school  of  fish.  In  this  model,  any  detection  not 
associated  with  the  single  real  target  will  be  called  a  false  alarm. 
Every  position  report  will  therefore  require  classification  as  either 
a  target  contact  or  a  false  alarm.  The  case  where  other  'trackable' 
objects  are  present  will  not  be  considered. 

False  alarms  will  be  considered  to  occur  according  to  a  space-time 
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Poisson  process  with  rate  n  [units  are  (length  x  time"  )].  This  re¬ 
sults  in  false  alarms  occurring  uniformly  over  an  area  with  time 
spacing  of  events  being  governed  by  a  Poisson  process.  This  type  of 
detection  process  can  be  assumed  to  arise  from  random  disturbances  to 
the  system  or  from  real  but  rare  detections  on  other  objects. 

The  assumption  will  also  be  made  that  there  is  no  discernable  dif¬ 
ference  between  false  and  real  detections.  This  assumption  can  be 
considered  highly  pessimistic  when  compared  to  any  existing  tracking 
system.  The  result  is  to  drive  the  average  tracking  time  toward  a 
lower  bound. 

We  have  used  four  parameters,  D,  n,  and  a,  to  characterize  the 
dynamics  of  the  tracking  problem.  Since  there  are  only  two  dimensions, 
time  and  length,  actually  involved,  we  can  summarize  the  problem  with 
two  dimensionless  summary  measures.  These  are: 

a  =  a/k'ZttiD  =  rate  ratio 
3  =  -  accuracy  ratio 


C.  TRACKER  MEMORY 


The  concept  of  memory  raises  several  possibilities  when  applied  to 
tracking  systems  in  general.  Qualitatively,  an  expanded  memory  would 
seem  desirable,  but  must  be  assumed  to  incur  an  additional  cost  to  the 
tracker.  Two  questions  addressed  in  this  report  are  how  much  memory 
is  needed  and  in  what  way  should  memory  be  used. 

First  we  address  the  simplest  tracking  problem  in  which  a  single 
decision  is  made  following  each  detection.  Tracking  begins  with  a 
known  real  target  detection.  Each  subsequent  position  report  is  evalu¬ 
ated  as  real  or  false  instantaneously.  If  evaluated  false,  the  detec¬ 
tion  will  be  forgotten  and  the  problem  will  continue.  If  evaluated 
real,  the  problem  will  be  updated  and  restarted  from  the  target's  new 
position  estimate.  Tracking  time  is  considered  to  stop  at  the  time  of 
the  tracker's  first  mistake;  i.e. 

1)  Tracker  accepts  a  false  alarm  as  real. 

2)  Tracker  rejects  a  real  detection. 

These  rules  apply  to  what  will  be  called  the  memoryless  tracker. 

Figure  1  shows  the  results  of  three  separate  tracking  sequences 
simulated  with  the  maximum  likelihood  memoryless  model  developed  in 
this  report.  The  time  of  the  tracker's  first  mistaken  identification 
is  marked  with  a  o  in  each  case.  Also  marked  with  a  o  is  the  last 
time,  within  the  allotted  simulation  time,  that  the  tracker  correctly 
identified  the  real  target.  From  the  sample  cases  shown  it  is  not 
clear  that  the  tracker  has  actually  lost  his  tracking  ability  after 


making  only  one  mistake.  The  tracker  demonstrates  an  apparent  ability 
to  reacquire  the  target  in  some  cases.  It  is  clear  though,  that  there 
is  a  point  in  time  beyond  which  error  will  generally  continue  to  in¬ 
crease.  Any  reduction  in  tracker  error  heyond  this  point  can  be  con¬ 
sidered  to  be  almost  exclusively  due  to  the  random  chance  that  the 
target  might  walk  back  into  the  tracker  rather  than  due  to  tracker 
efficiency.  A  problem  arises  however  because  of  the  fact  that  there 
is  some  positive  probability  that  at  some  time  in  the  future,  the 
tracker  might  accidentally  reacquire  the  target.  A  trivial  tracker 
could,  in  fact,  extend  the  tracking  time  indefinitely  by  accepting  all 
detections  that  occur  and  defining  tracking  time  to  end  at  the  time 
of  the  last  correct  identification.  For  this  reason,  we  will  use  E(T), 
the  average  time  till  the  tracker  makes  his  first  mistake,  as  our 
summary  measure  of  tracker  effectiveness  and  consider  it  a  useful 
lower  bound  on  actual  tracking  time. 

The  tracker  with  memory  will  face  the  same  task  and  its  performance 
will  be  judged  in  a  similar  manner.  Ue  look  first  at  the  case  where 
the  accuracy  ratio,  3,  equals  zero.  In  the  diffusing  target  scenario, 
future  target  positions  can  then  be  considered  to  be  independent  of  all 
but  the  real  target's  present  position.  Reports  received  prior  to  the 
most  recent  target  detection  will,  therefore,  not  be  retained  in 
memory.  Rather  than  demand  an  immediate  decision  on  subsequent  detec¬ 
tions,  however,  the  tracker  with  memory  will  be  given  the  opportunity 
to  file  a  finite  number  of  detections  in  memory  for  a  future  decision. 
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Following  each  detection  he  will  decide  whether  any  single  subset  of 
detections  on  file  constitutes  the  actual  target  track.  If  the  de¬ 
cision  is  yes,  memory  will  be  scrubbed  of  all  but  the  most  recent  posi¬ 
tion  estimate  and  the  process  will  begin  again.  When  no  track  has 
been  accepted,  the  tracker  will  retain  all  detections  on  file  until 
his  memory  capacity  is  exceeded.  As  additional  reports  are  received 
he  will  sequentially  scrub  detections  on  file  according  to  a  prede¬ 
termined  rule  to  be  discussed  later. 

In  the  situation  where  the  accuracy  ratio,  8,  is  not  equal  to  zero, 
the  problem  is  handled  in  essentially  the  same  way.  Successive  real 
target  detections  are  not  independent  of  previous  ones  however,  and 
will  require  special  handling.  This  case  will  be  dealt  with  in 
chapter  4. 

So  that  effectiveness  of  the  memoryless  tracker  and  the  tracker 
with  memory  may  be  compared,  tracking  time  will  be  defined  in  essen¬ 
tially  the  same  way  as  before.  The  tracker  with  memory  will  begin  at 
time  zero  with  a  known  real  target  detection  and  end  when  he  makes  his 
first  mistake.  Mistakes  occur  in  any  of  the  following  ways: 

1)  A  false  alarm  is  included  in  an  accepted  set  of  detections. 

2)  A  real  detection  is  omitted  from  an  accepted  set  of  detections. 
Tracking  time  is  stopped  at  the  time  associated  with  the  mistaken 
detection,  which  may  be  earlier  than  the  time  at  which  the  fatal 
decision  is  made. 
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In  the  next  chapter,  we  will  develop  the  theory  and  specific 
decision  rules  for  what  will  be  called  the  maximum  likelihood  tracker. 
To  simplify  the  analytic  solution,  measurement  error  will  not  be  con¬ 
sidered  in  that  chapter  or  the  next.  In  chapter  3  we  will  examine 
the  results  of  a  maximum  likelihood  tracker  simulation  and  compare 
the  results  with  those  of  Washburn  [3].  A  memory  capacity  will  be 
given  to  the  tracker  in  this  chapter  and  its  effect  will  be  analyzed. 
In  chapter  4  we  will  introduce  and  deal  with  measurement  error. 

Chapter  5  will  focus  on  the  conclusions  of  this  report. 


II.  MAXIMUM  LIKELIHOOD  DECISION  RULE 


In  classical  statistics  the  maximum  likelihood  estimator  is  typi¬ 
cally  used  to  estimate  the  value  of  a  distributional  parameter  that 
has  maximum  chance  of  producing  a  sample  result.  In  the  study  of  a 
time-sequenced  series  of  events  the  same  concept  can  be  applied.  The 
subject  of  this  chapter  will  be  to  develop  a  simplified  decision  rule 
that  reflects  the  concept  of  maximum  likelihood.  The  motivation  lies 
behind  the  idea  of  labeling  any  series  of  detections  in  such  a  way  as 
to  maximize  the  probability  of  that  set  of  labels  being  correct.  The 
subset  of  detections  labeled  real  in  a  series  of  detections  could 
also  be  called  the  target  track.  The  final  result  will  be  a  simple 
summary  measure  by  which  every  series  of  detection  labels  can  be  com¬ 
pared.  Only  the  most  likely  series  of  events  will  be  accepted. 

A.  SUMMARY  OF  TRACKING  RULES 

The  problem  of  evaluating  a  single  detection  is  particularly 
adaptable  to  a  maximum  likelihood  model .  We  need  only  compare  the 
probability  of  the  detection  being  real  against  the  probability  of 
its  being  false.  Given  the  characteristics  of  the  diffusing  target, 
this  problem  will  yield  easily  to  analysis. 

When  a  series  of  K  detections  must  be  considered,  the  problem 
escalates  in  difficulty  but  follows  the  same  basic  principles.  All 
possible  combinations  of  real  and  false  detection  labels  must  be 
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compared.  The  combination  of  labels  that  holds  the  greatest  chance  of 
being  correct  would  be  accepted  as  real.  The  label  combinations  would 
be  bounded  by  the  set  where  all  detections  were  labeled  false  as  well 
as  by  the  set  where  all  were  labeled  real.  All  combinations  within 
those  bounds  would  have  to  be  examined. 

With  this  type  of  tracker  it  would  be  theoretically  possible  but 
practically  unmanageable  to  retain  all  detections  on  file  for  future 
re-evaluation,  even  after  they  may  have  already  been  accepted.  For 
the  purposes  of  this  report,  the  tracker  will  not  be  permitted  to 
reconsider  past  decisions.  The  tracker  can  only  progress  from  the 
last  accepted  decision. 

The  magnitude  of  this  problem  increases  dramatically  as  additional 
detections  are  added  to  tracker  memory.  At  each  step  the  number  of 
possible  combinations  which  must  be  considered  actually  doubles. 
Accordingly,  the  size  of  tracker  memory  quickly  becomes  an  issue  of 
critical  concern. 

To  simplify  the  analysis  somewhat,  the  tracker's  circular  normal 
measurement  error  will  not  be  considered  in  this  chapter  or  the  next. 
The  inclusion  of  measurement  error  serves  to  complicate  computer  im¬ 
plementation  of  the  decision  process,  although  the  maximum  likelihood 
rules  remain  essentially  unchanged.  The  case  where  measurement  error 
is  not  zero  will  be  dealt  with  in  chapter  4. 

To  summarize  the  tracking  rule  that  will  be  developed  in  the 
remainder  of  this  chapter,  we  will  first  define  a  new  concept  for 
distance.  Any  two  detections  will  have  a  dimensionless  "distance" 
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between  them  equal  to  [R /(2Dt)]  +  1  n(. 2rrn  Dt/A)  where  R  is  the  physical 
distance  and  t  is  the  time  difference  between  detections.  A  "track" 
then,  will  simply  be  a  set  of  detections  connected  in  order  of  occur¬ 
rence,  and  the  most  likely  track  is  the  one  with  the  shortest  total 
distance.  Note  that  the  distance  can  be  negative  so  that  the  best 
track  is  not  necessarily  null.  In  fact,  the  null  solution  corresponds 
to  case  where  every  detection  is  evaluated  as  being  false.  The  rule 
presented  here  is  actually  a  special  case  of  the  tracking  scheme 
developed  by  Sittler  [2]. 

3.  DEVELOPMENT  OF  TRACKING  RULES 

To  begin  the  analysis  we  consider  a  discrete  approximation  to  the 
problem.  We  will  assume  time  is  discrete  with  intervals  of  length  a. 
Space  will  be  considered  to  be  made  of  discrete  sections  with  area  of 
size  5.  Figure  2  illustrates  a  possible  sequence  of  two  events  in 
space  and  time.  Although  figure  2  shows  only  one  dimension  of  space, 
x,  a  second  dimension,  y,  will  also  be  considered  in  the  mathematics 
that  follows.  With  two  detections  on  file  we  must  consider  labeling 
the  events  in  one  of  four  possible  ways 

1)  Detection  one  is  real,  two  is  false 

2)  Detection  one  Is  false,  two  is  real 

3)  Detection  one  is  real,  two  is  also  real 

4)  Detection  one  is  false,  two  is  also  false 

The  likelihood  of  any  of  these  label  patterns  or  tracks  can  be  ex¬ 
pressed  with  a  combination  of  the  following  expressions: 


Figure  2.  Sample  space-time  plot  for  two  detections. 


P[Real  detection  occurs  at  (X..  ,Y^ )  at  time  t^]  » 

(6/ ( 2ffDti )  expi-(Xf  2+Y . 2 )/[  2Dtf )  ] )  ( \A )  ( 1 ) 

PfFalse  alarm  occurs  at  (X^ ,Y^)  at  time  t^]  = 

(«S/A)(nM)  «  cnA  (2) 

The  first  factor  in  expression  (1)  is  simply  the  normal  density  func¬ 
tion  for  a  diffusing  target  multiplied  by  the  elementary  area  5.  It 
represents  the  probability  that  the  target  is  within  a  given  elementary 
area  at  a  specific  time.  The  variable  t.  is  actually  n^A,  where  is 
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the  number  of  time  intervals,  A,  since  the  last  detection  labeled  real 
in  the  series.  N  will  be  the  total  number  of  time  intervals.  A,  in¬ 
cluded  in  the  set.  The  word  'at'  is  being  used  loosely  to  indicate  a 
detection  within  an  elementary  space  and  time  interval.  and  Yi  are 
the  rectangular  distances  from  the  last  detection  considered  real  in 
the  series.  The  second  term  in  the  expression  is  the  Poisson  proba¬ 
bility  of  a  real  detection  occurring  in  a  specific  elementary  time  in¬ 
terval  A. 

Expression  {2}  is  the  false  alarm  equivalent  to  expression  (1). 

To  avoid  the  issue  of  boundary  conditions,  in  expression  (1)  we  must 
define  the  area,  A,  in  expression  (2),  to  be  large  enough  to  effective¬ 
ly  negate  any  possibility  of  the  target  reaching  a  boundary  before 
either  the  tracker  loses  track  or  is  able  to  reinitialize  the  problem. 

Referring  back  to  the  four  possible  label  combinations,  the  proba¬ 
bilities  for  each  of  these  events  can  be  expressed  as  follows: 

P-|  *  P[R.D.  in  <5.j  when  n^ -A=  t^]*P[F.D.  in  when  n2*^  =  t?! 

•  P[No  R.D .  in  (N-l  )A' s]*P[No  F.D.  in  (N-l)A's]  (3) 

?2  *  P[F.D.  in  o-j  when  n-j  *A  =  t-j]*P[R.D.  in  52  when  n2*A=  t 

•  P[No  R.D.  in  (N-l)A's]*P[No  F.D.  in  (N-l)A's]  (4) 

P-j  s  P[R.D.  in  5-|  when  n^*A  *  t^]*P[R.D.  in  62  when  n2*A*  1 2 J 

•  P[No  R.D.  in  (N-2)A's]*P[No  F.D.  in  (N)  A's]  (5) 
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P^  =  P[F.D.  In  51  when  n^ *A  *  t^]«P[F.D.  in  $2  when  n2*A  a  t2] 

•  P[No  R.D.  in  N  A's]  •  P[No  F.D.  in  (N-2)A's] 

Then  expressing  (3)  and  (5)  analytically  we  have: 

P1  =  C(XA)*6/(2ir0ti)exp[-(X12+Y12)/(2Dt1)]] 

•  [(nM)(<S/A)]Cl-AA]T/A_1  [1-nAA]™-1 


P3  -  [CXA)-5/(27rDt1)exp[-(X12+Y12)/(2Pt,)]] 

•  C(AA)(6/(2TrDt2)exp[-(x/:  Y.\ 

•  [1-XA]™‘2  [l-nAA]T/A 


where  T  *  N  A  is  the  total  tHs  being  considered.  The  remaining 
terms  can  be  expressed  in  similar  fashion. 

The  task  now  is  to  find  a  way  in  which  to  compare  the  various 
analytic  probabilities.  It  can  be  seen  immediately  that  each  of  the 
expressions  will  contain  a  number  of  equivalent  factors. 

Shifting  momentarily  to  the  general  case,  we  will  consider  the 
situation  where  there  are  K  detections  in  the  set,  m,  of  which  are 

a 

being  evaluated  as  real  in  track  (a)  and  mb  of  which  are  being  evalu¬ 
ated  as  real  in  track  (b).  The  likelihood  ratio  of  the  two  tracks  is 


[(AA)  a  TT  5/C2TDtat)  exp[-(Xa.2+Yai2]/C20taf]] 
i=l 


[CXA)"b  TT  5/C2iDtbj)  exp[-(Xbj2+Ybj2)/(2Dtbj)] 
=1  'J  *JJ 


.  C(S/A)K~ma  (nAA)K  (l-AA)™~ma  (l-nAA)™~(K~ma) 
CC6/A)K_mb  (nM)K"mb]  (l-AA)™'mb  (1  -nAA)T/A-CK-mb5 


(9) 


which  reduces  to 


m  a  K— m 

f C AA )  aJT*6/C27rDta.)  exp [-( Xfl i 2+Ya i 2 ) / C 2Dta ^ ) 3 C 6nA )  a 

.  1=1 _ 

Pb  mb 

m,  ,  K-m. 

[(AA)  b  Jj  5/(2Tr0tbj)  exp[-(XbjZ+Ybj/)/(2Dtbj)](5nA) 

j=l 


m, -m  m  -m. 

(1 -AA)  b  a  (1-nAA)  3  b 


(10) 


The  subscripts  on  X,  Y,  and  t  have  been  changed  to  indicate  that  there 
are  different  values  associated  with  the  terms  from  different  paths; 

X  .  is  the  X-coordinate  of  the  iTH  detection  on  track  (a),  etc.  Note 

a  I 

that  when  A  is  small  the  last  two  factors  in  equation  (10)  are  approxi¬ 
mately  equal  to  the  value  one.  We  will  make  this  approximation,  and 
by  further  cancellation  we  obtain 
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(11) 


(A/nJ^TT  exp[-{Xa.2+Ya1.2)/(2Dta.)] 

_ M _ 

m  mb 

(X/n)  bjy  1/(2irDtbj)  exp[-(Xbj2+Ybj2)/(2Dtbj)] 
j=1 


which  suggests  a  useful  decision  rule.  Note  that  the  quantities  A, 

A,  and  o  do  not  appear.  Each  term  in  the  expression  is  either  known 
or  can  be  measured.  A  decision  rule  would  be  to  pick  track  (a)  if  the 
ratio  was  greater  than  one,  or  pick  track  (b)  if  the  ratio  was  less 
than  one.  The  comparison  could  be  carried  out  for  all  possible 
tracks  in  the  set  of  detections  with  the  winner  being  declared  most 
likely.  The  method  is  intuitively  reasonable  and  easily  implemented 
by  computer  as  long  as  the  maximum  number  of  possible  tracks  remains 
small . 

For  simplification,  it  is  convenient  to  use  the  logarithmic  form 
of  equation  (11).  Let  Li=ln  Pi.  Then  Pi >P j  if  and  only  if  Li <L j . 

The  rule  then  would  be  to  pick  the  track  with  the  smallest  value  Li. 
The  final  summary  measure  for  track  (a)  will  be 

m 

a 

La  “2  ClXalHi2,/<2Dtai1  *  1  "(2™»VX)]  (12) 

i*l 

The  criterion  for  track  acceptance  will  be,  accept  the  track  with  the 
smallest  value  Li.  The  case  where  all  detections  are  labeled  false 
will  have  the  value  zero.  In  this  way  we  are  actually  choosing  to 


label  the  observed  series  of  events  with  the  pattern  of  detection  types 
that  is  most  likely  to  have  occurred. 

Note  that  the  first  term  in  equation  12  is  always  positive.  The 
second  term  will  be  negative  only  when  the  ratio  ( 2im  DTi  A )  is  less 
than  one.  The  implication  is  that  there  is  a  limit  in  time  beyond 
which  a  single  detection  can  not  be  accepted.  With  no  memory,  the 
tracker  will  stall  if  he  has  not  accepted  a  detection  by  a  given  time 

limit.  With  a  memory  capacity,  the  tracker  is  able  to  weight  the 

various  combinations  of  time  and  space  differences  between  detections. 
In  this  way  he  may  still  accumulate  an  acceptable  track. 

As  an  example  we  can  look  back  at  the  situation  developed  from 
figure  2.  The  four  sumary  measures  would  be 

1)  L,  *  (X112+Y112)/(20t11)  +  ln(2im0t11A)  (13) 

2)  L2  »  (X212+Y212)/(2Dt21)  +  ln(2^Dt21A)  (14) 

L3  s  ^X31  +Y31  ) / C 2Dt3-j )  +  lr>(X32  +Y32  )/(2Dt32) 

+  ln(2miDt3i  A)  +  1  n(2-nnDt32A)  (15) 

4)  L4  =  0  (16) 

The  trackers  rule  would  be  to  accept  the  track  with  the  smallest  value, 
update  the  memory  file  accordingly,  then  wait  for  the  next  detection 
to  occur. 


23 


III.  TRACKING  WITH  NO  MEASUREMENT  ERROR 


For  the  purpose  of  examining  the  value  of  a  new  tracking  model,  it 
is  useful  to  make  comparisons  with  a  model  already  in  existence.  In 
this  section  we  will  examine  results  of  simulating  the  maximum  likeli¬ 
hood  tracker  and  compare  the  results  with  those  of  Washburn  [3],  who 
has  developed  an  optimized  tracking  rule  for  the  memoryless  tracker 
with  no  measurement  error.  We  will  then  introduce  memory  to  the  maxi¬ 
mum  likelihood  model  and  note  the  improved  results. 

A.  THE  MEMORYLESS  TRACKER 

Examining  the  same  problem  presented  here,  Washburn  [3]  arrives  at 
an  analytic  solution  which  produces  a  tracking  rule  based  on  detection 
range  as  a  function  of  time.  Given  the  dynamic  parameters  of  the 
problem  and  a  detection  time,  his  rule  will  generate  a  range  from  the 
last  target  position  inside  which  the  new  detection  should  be  accepted 
as  real.  Should  the  detection  fall  outside  that  range  circle  it  would 
be  rejected  as  false.  The  method  is  relatively  easy  to  implement 
requiring  only  a  single  function  evaluation  for  each  detection  that 
occurs.  Based  on  this  algorithm  he  is  able  to  identify  a  dimensionless 
curve  for  XT,  (detection  rate  x  average  tracking  time),  as  a  function 
of  the  rate  ratio  a.  This  curve  is  reproduced  in  figure  3.  He  has 
shown  that  these  are  maximum  expected  tracking  times  for  the  memoryless 
tracker  with  no  measurement  error. 
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The  continuous  time  Fortran  simulation  listed  at  the  end  of  this 
report  was  designed  to  examine  the  effectiveness  of  the  maximum  likeli¬ 
hood  decision  rule.  The  dynamic  parameters  of  the  program  include  the 
component  parameters  of  the  rate  ratio,  2,  and  the  accuracy  ratio,  3, 
as  well  as  a  memory  input  which  allows  for  tracker  memory  set  between 
one  and  ten  detections.  The  case  where  the  tracker's  maximum  memory 
capacity  is  equal  to  one  is  the  equivalent  of  Washburn’s  memoryless 
tracker.  Using  the  maximum  likelihood  memoryless  tracker  model,  it  is 
possible  to  generate  a  graph  similar  to  that  of  Washburn's.  Point 
estimates,  based  on  1000  iterations  at  each  a  level  are  marked  with  a 
q  in  figure  3.  Table  I  shows  Washburn’s  XT,  the  maximum  likelihood 
point  estimates,  a  95%  confidence  interval  for  the  point  estimates,  and 
the  sample  standard  deviation  for  the  maximum  likelihood  results. 

The  results  indicate  that  tracking  time  for  the  maximum  likelihood 
memoryless  tracker  is  roughly  exponential  in  nature,  with  the  sample 
mean  approximately  equal  to  the  sample  standard  deviation.  More  im¬ 
portantly,  the  times  are  nearly  as  good  as  those  of  Washburn.  Wash¬ 
burn's  model  was  specifically  designed  to  maximize  tracking  time  as 
it  is  defined  in  his  report.  Using  the  same  definition  for  tracking 
time,  we  are  able  to  nearly  duplicate  his  results  with  the  simplest 
version  of  the  maximum  likelihood  tracker. 

B.  THE  TRACKER  WITH  MEMORY 

Having  established  a  tracking  algorithm  that  is  nearly  the  equal  of  an 
optimized  memoryless  tracker,  our  goal  will  now  be  to  improve  upon  the 
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TABLE  I.  Dimensionless  tracking  value,  AT,  as  a  function  of 
rate  ratio,  a;  3=0. 


a 

Washburns 

AT 

Maximum 
Likel ihood 
AT 

95%  Confidence  Interval 
Lower  Upper 

Bound  Bound 

Sampl  e 
Standard 
Deviation 

1 

1  .4 

1  .21 

1  .15 

1  .27 

1 .02 

2 

2.35 

2.33 

2.22 

2.44 

1 .83 

3 

3.8 

3.52 

3.34 

3.70 

2.94 

4 

5.5 

5.22 

4.93 

5.51 

4.65 

5 

7.35 

6.95 

6.56 

7.34 

6.22 

6 

9.45 

9.35 

8.81 

9.89 

8.72 

7 

11.9 

11.28 

10.62 

11  .94 

10.60 

8 

14.4 

14.04 

13.12 

14.95 

13.36 

9 

17.1 

17.20 

16.11 

18.29 

17.52 

10 

20.1 

20.08 

18.92 

21  .24 

18.69 

TABLE  II. 

Dimensionless  tracking  value,  AT,  as  a 
memory;  a =4,  3=0. 

function  of 

Memory 

AT 

95%  Confidence 
Lower 

Bound 

Interval  on  AT 
Upper 

Bound 

%  Improvement 

100,(A5'.39~~) 

1 

5.39 

5.08 

5.70 

— 

Washburn's  AT  =  5.45 

2 

5.53 

5.23 

5.83 

2.3 

3 

5.84 

5.54 

6.14 

8.3 

4 

5.86 

5.54 

6.18 

8.7 

5 

5.82 

5.49 

6.15 

7.9 

6 

5.90 

5.57 

6.23 

9.5 

7 

5.90 

5.59 

6.21 

9.5 

8 

5.58 

5,28 

5.88 

3.5 

9 

6.13 

5.78 

6.48 

13.7 

10 

5.90 

5.58 

6.22 

9.5 
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model  fay  introducing  a  memory  capacity.  As  stated  in  part  A  of  this 
chapter,  the  simulation  used  in  this  report  was  run  with  as  many  as 
ten  detections  on  file.  At  each  detection  the  model  evaluates  every 
possible  label  combination  for  the  detections  in  memory  and  accepts 
the  most  likely.  Note  that  the  track  selected  will  either  include 
the  most  recent  detection  or  be  null.  In  accordance  with  the  tracking 
rules,  when  a  non-null  track  is  accepted  all  but  the  most  recent  posi¬ 
tion  estimate  are  scrubbed  from  memory.  The  result  is  a  type  of  time 
sequenced  problem  where  the  number  of  label  combinations  to  be  examined 
actually  doubles  as  each  new  detection  is  added  to  memory.  For  the 
purpose  of  this  report,  a  memory  capacity  of  ten  detections  was  found 
to  sufficiently  illustrate  the  benefit  of  memory.  Further  increases 
in  memory  capacity  are  made  at  a  substantial  cost  in  computer  process¬ 
ing  time.  As  will  be  shown  later  in  this  report,  the  small  increase 
in  tracking  time  is  not  likely  to  be  worth  the  extra  effort  required 
for  more  memory.  There  appears  to  be  an  upper  bound  on  tracking 
ability.  The  bound  is  nearly  reached  with  a  memory  capacity  less  than 
ten.  As  the  capacity  increases  tracking  time  asymptotically 
approaches  the  upper  bound. 

Because  of  the  limited  storage  capacity,  the  possibility  does  exist 
that  memory  capacity  might  be  exceeded  before  a  detection,  or  set  of 
detections,  is  accepted.  In  this  case  we  will  have  to  introduce  a 
procedure  whereby  detections  can  be  selectively  purged  from  memory. 

The  question  of  where  to  begin  memory  cleansing  stimulates  several 
possible  approaches.  Not  only  individual  detections,  but  combinations 
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of  several  detections  could  be  considered  to  determine  which  might 
offer  the  greatest  chance  of  future  acceptance.  For  simplicity,  in 
this  report  the  detections  will  be  considered  individually,  and  the 
single  most  unlikely  detection  will  be  scrubbed. 

Looking  now  at  a  specific  example,  we  begin  to  observe  the  effect 
of  memory.  In  the  case  where  a=4,  3=0,  the  incremental  effects  of 
memory  are  illustrated  in  figure  4.  A  curve,  based  on  1000  iterations 
at  each  memory  level,  has  been  fitted  between  the  points  to  roughly 
approximate  the  shape  of  the  empirical  solution.  Although  not  smooth, 
the  curve  does  indicate  a  significant  gain  in  tracking  time  as  memory 
capacity  is  initially  expanded  to  several  detections.  The  curve  then 
seems  to  flatten  somewhat,  indicating  a  diminishing  return  for  memory. 
This  flattening  of  the  curve,  found  in  every  sample  run,  provided 
justification  for  limiting  memory  capacity  to  ten  detections.  In 
fact,  if  computer  access  were  limited,  a  memory  capacity  of  three  or 
four  detections  might  prove  sufficient.  Table  II  lists  the  point 
estimates,  95%  confidence  intervals,  and  the  per  cent  improvement  over 
base  case,  memory  =  1  . 

Although  we  have  reported  only  one  case,  the  effects  of  memory 
were  found  to  be  universal.  Various  runs  with  a  ranging  from  one  to 
ten  showed  the  maximum  likelihood  tracker  could  consistently  outperform 
Washburn's  memoryless  tracker  when  equipped  with  a  memory  capacity  as 
small  as  two  or  three  detections.  In  the  next  chapter  we  will  examine 
the  method  by  which  the  maximum  likelihood  tracker  is  able  to  handle 
the  problem  of  measurement  error. 


IV.  TRACKING  WITH  MEASUREMENT  ERROR 


Measurement  error  in  the  detection  system  is  a  common  problem  for 
trackers.  In  this  report  the  error  is  assumed  to  be  circular  normal. 
This  distribution  is  commonly  used  in  modeling  both  because  of  its 
appropriateness  for  many  existing  systems  and  because  of  its  mathemati¬ 
cal  tractabil ity.  In  this  section  we  will  discuss  the  effects  of 
measurement  error  on  the  maximum  likelihood  tracker  and  expand  the 
model  to  deal  with  it  accordingly.  In  addition,  we  will  examine  the 
results  of  a  computer  simulation  in  which  measurement  error  was 
included . 

To  illustrate  the  added  dimension  that  measurement  error  imposes, 
we  can  first  look  at  figure  5.  The  probabilistic  term  for  real  target 
position  reports  can  no  longer  be  defined  with  the  simple  bivariate 
normal  motion  equation  used  in  chapter  2.  Instead,  we  will  define 
a  new  expression  which  accounts  for  both  the  motion  and  the  error  be¬ 
tween  subsequent  position  reports.  In  figure  5,  R.  i+1  will  be  the 
vector  distance  between  the  target's  actual  position  at  times  i  and 
i+1 .  R.  will  be  the  measurement  error  associated  with  the  i TH  detec¬ 
tion.  The  vector  distance  between  real  target  detections  will  be 
called  Q..  In  this  case  we  have 

Q1  =  “  R1  +  R1 ,2  +  R2 

Q2  =  -  R2  +  R2,3  +  R3 

Qi  “  ‘  Ri  +  Ri,i+1  +  Ri+1 
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Figure  5.  Sample  set  of  target  position  reports  and 
associated  target  positions. 


Since  the  problem  is  symmetric  in  the  X  and  Y  dimension,  for  sim¬ 
plicity  we  can  examine  a  one  dimensional  version  of  the  same  problem. 
Since  the  X  and  Y  components  of  target  motion  and  error  are  independent 
and  parameterized  identically,  we  will  be  able  to  combine  terms  to  form 
a  single  solution  based  on  the  identically  distributed  but  independent 
resul ts . 

In  this  case  the  one  dimensional,  horizontal  distance  between  de¬ 
tections  will  be  called  X^ .  The  error  associated  with  the  iTH  detec¬ 
tion  will  be  called  Ei ,  and  the  distance  between  actual  target  position 
i  and  i+1  will  be  V.. .  A  series  of  measurements  can  be  represented  as 
follows : 
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X1 =  -  E1 +  V1 +  h 


X2  =  "  E2  +  V2  +  E3 


Xi  ’  -  Ei  +  Vi  +  Ei+1 


We  know  the  E  and  V  variables  are  all  independent  with  means  equal  to 

2 

zero,  that  the  variance  of  is  equal  to  a  ,  and  that  the  variance 
of  V.j  is  equal  to  Dt.. .  However,  the  individual  measurements,  X.,  are 
not  independent.  Since  E2  is  common  to  both  X^  and  Xj,  E [X-j  X3J  = 

-  a2,  and  more  generally  E[X^  X.+^3  =  -  Since  E[X2]  =  2cr2+Dt.]. , 
the  covariance  matrix  of  U3  =  CX-j.X-.Xj)^,  for  instance,  will  be 


-a 

0 

2a2+  Dtz 

-a2 

2 

-a 

2o2+  Dt 

This  matrix  retains  its  form,  expanding  only  along  the  diagonals,  as 
the  dimension  is  increased.  Since  the  vector  IL  =  (X^,  X2 , . . - ,  Xi)t 
is  multivariate  normal,  the  probability  density  function  of  1).  will  be 

f(U.)  =  y'lZ^ ,  _1)  exp[-  y  J.U] 
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where  and  |  |  denotes  determinant  [4].  3y  symmetry, 

the  same  probability  density  function  holds  for  the  vector  of  incre¬ 
mental  measurements  in  the  vertical  direction,  Wi  =  (Y^ ,  Y2>...,  Y^)t. 
Since  U.  and  are  independent,  the  joint  probability  density  function 
is  obtained  by  multiplication: 

f(U.,W.)  »  (M1  |E1  i  )_1  exp[-J./2] 

where  J.W  =  E*"1  and  Ji  =  J.U  +  J^.  This  then  will  replace 

the  term  used  in  the  decision  rule  developed  in  chapter  2  .  In  final 
form,  the  new  expression  will  be: 

Li  -  \  Ji  +  InCUimA)1  |Ei',] 

Note  that  this  is  the  equivalent  of  (12)  from  chapter  (2)  when  a 
equals  zero. 

Although  the  decision  rule  looks  relatively  simple  when  expressed 
in  matrix  notation,  computer  implementation  of  the  rule  creates 
another  problem.  To  avoid  the  burden  of  recomputing  E”1  and  jzj  for 
each  of  the  tracks  examined,  it  is  necessary  to  develop  an  iterative 
procedure  by  which  a  summary  term  can  be  identified  for  each  of  the 
tracks  examined.  Using  a  type  of  Kalman  filter,  developed  mathe¬ 
matically  in  Appendix  A,  it  is  possible  to  obtain  a  position  estimate, 
track  value  CL. ^ ) ,  and  gain  matrix  for  each  of  the  tracks  examined.  If 
none  of  the  tracks  are  subsequently  accepted,  the  terms  can  be  saved 
and  used  when  evaluating  the  next  detection  that  occurs.  The  procedure 
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is  completely  iterative  and  reduces  to  the  procedure  developed  in 
chapter  (2)  when  3  equals  zero.  A  flow  chart  for  the  algorithm  is 
outlined  below  with  S  =  Dti  and  d  =  \z\ .  See  Appendix  A  for  the 
definition  of  v,  x,  and  y. 

J  =  0 

y  a  co 

x  =  0 
y  =  0 

d  =  1 

i - >  Observe  X,  Y,  and  S2 

i  ?  ?  4. 

|  v  =  2a  +  S  -  a  /v 

i 

d  =  d  •  v 

J  =  J  +  [(x-X)2  +  Cy«Y)2]/v 

j  x  »  a2(x-X)/v 

*  ? 

|  y  =  cT(y-Y)/v 


To  test  the  effectiveness  of  the  revised  decision  rule,  we  can  use 
the  Fortran  simulation  listed  at  the  end  of  this  report.  For  a  sample 
case  where  a=4  and  3=1,  the  sample  results  from  1000  iterations  at 
each  memory  level  are  graphed  in  figure  6.  Point  estimates,  95*  con¬ 
fidence  intervals,  sample  standard  deviation,  and  the  per  cent  improve¬ 
ment  over  base  case  (memory  =  1)  are  listed  in  Table  III.  The  results 
are  as  we  might  intuitively  expect.  Memory  capacity  again  appears  to 
benefit  the  tracker,  but  the  gain  appears  to  level  off  as  memory  is 


* 
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TABLE  III.  Dimensionless  tracking  time,  AT,  as  a  function  of 
memory;  a=4 ,  3=1  . 

95"  Confidence  Interval  Sample 


Memory 

At 

Lower 

Bound 

Upper 

Bound 

Standard 

Deviation 

Improve¬ 

ment 

1 

2.18 

2.05 

2.31 

2.04 

— 

2 

2.26 

2.13 

2.39 

2.13 

3.7 

3 

2.36 

2.22 

2.50 

2.24 

8.3 

4 

2.44 

2.29 

2.59 

2.37 

11  .9 

5 

2.55 

2.39 

2.71 

2.55 

17.0 

6 

2.47 

2.32 

2.62 

2.39 

13.3 

7 

2.47 

2.32 

2.62 

2.35 

13.3 

8 

2.39 

2.24 

2.54 

2.34 

9.6 

9 

2.60 

2.45 

2.75 

2.37 

19.3 

10 

2.51 

2.37 

2.65 

2.30 

15.1 

TABLE  IV.  Dimensionless  tracking  time,  AT, 
memory;  a=8;  3=1 . 

as  a  function 

of 

Memory  AT 

95"  Confidence 
Lower 

Bound 

Interval 

Upper 

Bound 

Sample 

Standard 

Deviation 

O! 

•0 

Improve¬ 

ment 

1  6.17 

6.05 

6.29 

6.04 

— 

2  6.53 

6.40 

6.66 

6.46 

5.8 

3  6.62 

6.49 

6.75 

6.55 

7.3 

4  6.58 

6.45 

6.71 

6.39 

6.6 
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increased  beyond  four  or  five  detections.  Compared  to  the  case  where 
8=0,  the  per  cent  improvement  in  tracking  time  appears  to  be  better 
although  the  upper  bound  on  tracking  time  is  much  less. 

Due  to  the  large  sample  standard  deviation,  the  95X  confidence 
intervals  are  also  large.  Consequently,  even  with  1000  iterations, 
the  point  estimates  do  not  produce  the  smooth  curve  we  might  have 
hoped  to  find  for  figure  6.  To  more  closely  approximate  the  nature 
of  the  curve,  we  will  examine  the  results  of  10,000  iterations  of  a 
test  case  where  3=8  and  3=1.  Memory  will  range  from  (1)  to  (4) 
detections.  The  results  are  graphed  in  figure  7,  and  listed  in 
table  IV.  As  expected,  the  results  begin  to  outline  the  shape  of  a 
smooth  curve,  rising  steeply  at  first,  then  tapering  off  to  some 
constant  value. 
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V.  SUMMARY 


First,  looking  at  the  case  where  there  is  no  measurement  error 
(3=0),  we  have  seen  that  the  maximum  likelihood  model  presented  here 
is  actually  a  special  case  of  the  maximum  likelihood  model  developed 
by  Sittler  [2].  In  this  model,  when  we  examine  a  set  of  detections,  a 
dimensionless  "distance"  is  assigned  to  every  possible  combination  of 
real  detections.  The  distance  is  based  on  a  maximum  likelihood  solution, 
and  the  combination  of  detections  with  the  shortest  total  distance  is 
accepted  as  the  target's  real  track.  The  method  is  intuitively 
reasonable  and  relatively  easy  to  implement  with  a  computer. 

When  this  method  is  compared  to  the  optimal  memoryless  tracker 
developed  by  Washburn  £3],  we  find  that  the  maximum  likelihood  memory¬ 
less  tracker  is  nearly  the  equivalent.  When  furnished  with  even  a 
limited  memory  capacity,  we  find  that  the  maximum  likelihood  tracker 
quickly  surpasses  the  optimal  memoryless  tracker  in  ability.  However, 
the  return  of  increased  tracking  time  for  additional  memory  diminishes 
quickly.  There  appears  to  be  an  upper  bound  on  tracking  f'me  which  is 
nearly  reached  with  a  memory  capacity  of  four  or  five  detections.  In 
all  cases  the  distribution  of  tracking  times  was  found  to  be  roughly 
exponential  with  the  mean  approximately  equal  to  the  standard  deviation. 

In  the  case  where  measurement  error  was  not  equal  to  zero,  we  were 
required  to  first  develop  a  type  of  Kalman  filter  to  handle  the  process 
of  obtaining  an  iterative  expression  for  the  "distance"  between 


detections  in  a  set.  Once  developed,  we  saw  that  the  case  where  3*0  was 
really  a  special  case  of  this  more  general  problem.  The  tracking  times 
still  appear  to  follow  an  exponential  distribution,  and  the  payoff  for 
memory  seemed  to  increase  as  3  approached  a. 

The  question  "How  much  memory  is  enough?",  is  particularly  relevant 
to  the  available  computational  power  as  well  as  the  importance  of  minor 
improvements  in  tracking  times.  It  appears  as  though  we  can  approach 
the  bound  with  as  few  as  four  or  five  detections;  however  some  small 
gain  in  tracking  time  will  still  be  realized  beyond  that  memory  level. 
The  cost  of  obtaining  the  small  increase  in  tracking  time  is  a  marked 
increase  in  computer  processing  time. 


i 
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APPENDIX  A 


DEVELOPMENT  OF  KALMAN  FILTER 


Me  will  consider  a  one  dimensional  system  in  which  position  measure¬ 
ments  are  taken  sequentially.  Target  motion  between  measurements  is 
assumed  to  be  the  result  of  normal  diffusion  and  will  be  represented  by 
V^,  where  ~  N(0,  SN  }  and  =  Dt^.  Each  position  measurement  will 

also  have  an  associated  measurement  error,  EN,  which  will  also  have  a 

2 

normal  distribution;  i.e.  E^  -  N(0,  c  ).  The  distance  between  succes¬ 
sive  measurements  will  be; 


Z„  =  E. 


VN  * 


-N-l 


N  =  1,  2, 


The  symbol  has  replaced  XN  of  the  text  (CH-4)  because  we  will  have  a 
special  use  for  the  symbol  X^. 

Let  be  the  covariance  matrix  of  =  [Z^ ,  ....  and 

Un .  Me  wish  to  find  an  iterative  means  of  producing  . 

To  begin  let  X^  =  [E^ ,  E^]*,  so  that  XN+]  =  5XN  +  WN;  where 

*  =  [l  o]  ’  MN  =  ^EN+1’  *  and  C0V^WN^  =  [o  o]  ■ 

Me  also  have  Z^  =  HX^  +  V^;  where  H  =  [1,  -1].  Now  let  the  NxN 
matrix  L^  be  a  whitening  filter;  i.e.  L^  is  white  noise.  Then 

!NxN  =  C0V£LN  un3  =  ln  en  ln  . 

Using  matrix  manipulations  we  obtain 

L  (L  ^  s  <~ 

LN  UN  1  “N 

and 
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therefore 


UNt^LNt  LN*  UN  =  UN  “  UN  =  JN  * 

Now  let 

LN  UN  =  Q2’  •••>  so  that 

N 

JN  =  ^LN  V*  ^LN  V  =  S  Qi  • 

i  =1 

A1  so  1  ot  2^  -  E[Z^  j  Z«|  ,  •  • « i 

VN  '  var  [ZnIZ1-  ••'•  ZN-1J  • 

Then 

Zn  ”  ZN 

— — —  is  independent  of  Z^ ,  . . . ,  ZN  ^  ,  normal,  and 

/VN 

has  unit  variance;  i.e. 

.  .  ZN  '  ZN 

Qn  "  — 

K  VN 

At  this  point,  our  intention  will  be  to  use  Kalman  filtering  to 

N  2 

obtain  an  iterative  expression  for  z^  and  v^.  Since  =  Z  and 

can  be  measured,  when  z^  and  v^  are  available  we  have  an  iterative 
expression  for  .  The  remainder  of  this  appendix  will  closely  follow 
the  Kalman  filtering  technique  used  by  Gel b  [1],  So  that  our  notation 
will  conform  to  his,  we  will  first  define  the  following  terms; 

XN(-)  3  eCxn  zi  •  •  ••  •  zn-i1  =  [_zj 
xN(+)  =  E[xN;zlf  ...,  ZN] 
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P^(-)  s  cov[X^j  Zi ,  .... 
P^(+)  =  *-0 v C Xj^ i Z-^  s  •  •  •  >  Z^J 


2  2  2 
To  begin,  let  =  var  [E^  i ^  ,  ....  ZN_.|]  >  so  that  o^  =  a  and 


V->  - 


o'  0 

Vj 


.  Then 


VN  »  H  PN(-)  Ht  +  SN2  =  a2  +  aN2  +  SN2 


The  Kalman  gain  matrix,  1*^,  is  given  by  K^  * 
the  Kalman  format,  we  can  update  the  current 
mate,  xNC+) ,  as  follows; 


a  ^VN  .  Following 


■o2/v 


error  esti- 


xN(+)  =  y-)  ♦  knczn  -  h  iNc-)] 


0 

^ZN  J 


‘2/vn1 

.-o2/v.J 


[V°"ZN] 


(0!/*n1Un-zn) 

.-ZN  -Ca2/ViN),CZN-zN)_ 


Since  x^.j(-)  =  iix^(+)  and  z^+-j  -  H  x^+-j(-) ,  z^^  =  H  i  x^(+)  - 
-ta2/vN) (.Z^-z^)  which  gives  the  desired  iterative  solution  for  z^.  We 
continue  working  toward  a  solution  for  v^  by  first  updating  the  co- 
variance  matrix; 


•  (I-KH)  PN(-) 


C I -KH ) 


where  qN 


2 


2  4. 

*  o  -  a  /v 

2 2,  \ 

=  a  (l-o  /vN) 
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Regardless  of  b,  c,  and  d,  we  can  compute  P^C-l  as  follows; 


PN+1(-)  s  $  PNC+)$  +  cov[WN] 


so  that  a 


2 

N+1 


At  this  point  we  have  completed  the  one-dimensional  analysis.  The 
computations  are  in  the  following  order  (an  intermediate  computation 
of  aN  has  been  eliminated!: 

Setting  v0  =  ”  to  initialize  the  problem 
we  have;  =  0 

V1  '  ♦  s/+  ®HW/v0) 

0,  ■  Z,2/v, 

z2  *  a2(z1  -  Z1 J/v1 

v2  =  a2  +  S22  +  a2(l-o2/v1) 

^2  =  Oi  *  i^2~zZ^  ^v2 


By  using  a  flow  chart  we  can  eliminate  the  variable  subscripts. 
The  procedure  will  be  as  follows; 

v  a  00 

z  =  0 
J  «  0 
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->  observe  Z  and  S 


o  2 
v  -  2a"  + 


4. 

a  /v 


i  j  =  j  +  (z  -  zr/v 
!  ? 

[j  *  a  (z  -  Z)/v 

When  we  examine  a  two  dimensional  ,  symmetric  version  of  the  same 
problem,  as  considered  in  the  text,  we  are  actually  looking  at  two 
independent  but  identical  problems.  In  this  case  we  solve  for 
JNU  =  U.,1  Z  _1  "  — '  ’  V  -  "  1  -  "1 


N  UN  ,  and  =  VN  Z^  using  the  procedure  developed 


in  this  appendix.  Because  of  the  symmetric  parameters,  we  can  combine 

terms  to  form  a  single  value  J.  In  a  flow  diagram  we  can  compute  J  as 

2  2 

follows;  J  =  J  +  [(x-X)  +  (y-Y)  ]/v  where  J  is  initially  set  equal  to 

zero.  This  term  can  then  be  used  in  a  bivariate  normal  probability 
distribution  used  to  represent  successive  real  detection  positions. 

The  entire  flow  chart  for  the  two  dimensional  case  is  shown  in 
chapter  (4)  of  the  text. 

Using  this  formulation,  we  can  also  find  a  way  to  iteratively 
calculate  the  value  of  jz.|.  As  previously  defined  QN  =  LNUN  and 


cov[Q]  =  E[Qn  Q^]  =  lM  ZN  »  INy^,so  that  cov[UM] 


’N  N  UN 


NxN  1 


'N1 
L  t)~1 

ln  un  ' 


•  th 


Since  will  be  a  lower  rectangular  matrix  with  the  i  diagonal 
element  equal  to  l//v7,  we  have; 


iave;  nr 

,LNi  o/y-rrv.  =  ;l 


t. 

N  1 


i=l 


and  therefore; 


;LN  1 


N 

TT 

i=l 


-  5  ^  vi 


since 
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;AB|  *  ■ A S  |B|  and  |A’^j  *  ] A j “1 .  This  then  will  be  the  iterative 
solution  for  |ZN|,  as  is  also  shown  in  the  flow  diagram  of  chapter  4. 
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